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On the Exact and Approximate Eigenvalue 
Distribution for Sum of Wishart Matrices 

Santosh Kumar, Gabriel Fernando Pivaro, Gustavo Fraidenraich, and Claudio Ferreira Dias 


Abstract 

The sum of Wishart matrices has an important role in multiuser communication employing multiantenna elements, such as 
multiple-input multiple-output (MIMO) multiple access channel (MAC), MIMO Relay channel, and other multiuser channels 
where the mathematical model is best described using random matrices. 

In this paper, the distribution of linear combination of complex Wishart distributed matrices has been studied. We present a 
new closed form expression for the marginal distribution of the eigenvalues of a weighted sum of K complex central Wishart 
matrices having covariance matrices proportional to the identity matrix. The expression is general and allows for any set of linear 
coefficients. 

As an application example, we have used the marginal distribution expression to obtain the ergodic sum-rate capacity for the 
MIMO-MAC network, and the cut-set upper bound for the MIMO-Relay case, both as closed form expressions. 

We also present a very simple expression to approximate the sum of Wishart matrices by one equivalent Wishart matrix. 

All of our results are validated by means of Monte Carlo simulations. As expected, the agreement between the exact eigenvalue 
distribution and simulations is perfect, whereas for the approximate solution the difference is indistinguishable. 

Index Terms 

Sum of Wishart matrices, eigenvalue distribution, multiple-input multiple-output, ergodic sum capacity, Meijer-G function. 

I. Introduction 

A. Random matrices and MIMO single-user relation 

R andom matrix theory has evolved into a truly multidisciplinary subject with its applications in helds as varied as commu¬ 
nication theory, quantum transport, quantum chromodynamics, quantum information theory, string theory, econophysics, 
number theory, etc. m. It is possible to represent the operators relevant to study physical systems in matrix form and use its 
properties to tackle difficult problems. Communication theory is one of the prominent areas on which random matrix theory 
has had huge impact. Random matrices gained special attention in wireless communication field after the works of Winters 
0, Foschini El, and Telatar 0. They have shown that the use of multiple antennas could enhance capacity in systems with 
limited bandwidth. In all cases, mathematical tools for matrices were employed for the analysis. 

In Telatar’s paper 0, the multiple-input multiple-output (MIMO) point-to-point ergodic channel capacity has been found. 
He has shown that instead of dealing with the joint probability density function of a Wishart distributed matrix 0, which is 
not easy to handle even for small dimensions, it suffices to use the joint eigenvalue probability density function given by James 
0. Such a simplihcation is possible in view of the unitarily-invariant nature of the channel capacity and other metrics which 
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2 


are usually employed to characterize the MIMO systems. Telatar’s pioneering work was the first one to establish a connection 
between MIMO communication and random matrix theory. Since then, there has been a great deal of interest in exploring and 
comprehending the properties of Wishart matrices. 

As a special case of Hermitian matrices, Wishart matrices arise in scenarios where MIMO systems are subject to Rician 
or Rayleigh fading. As indicated above, the performance of MIMO systems can be statistically predicted with the aid of 
eigenvalues distribution of Wishart matrices a, Q. For example, the channel matrix in a MIMO system relates to a Wishart 
matrix, whose eigenvalue statistics then leads to the knowledge of the ergodic capacity of the MIMO channel 0. On the other 
hand, the distribution of the largest and smallest eigenvalue can be used to analyze the performance of MIMO maximal ratio 
combining systems and MIMO antenna selection techniques, respectively ITJ. In JS), the authors have shown that the symbol 
error rate (SER) performance of MIMO systems employing multichannel beamforming in arbitrary-rank Ricean channels is 
dominated by the subchannel SER corresponding to the minimum channel singular value. Their results are based on marginal 
ordered eigenvalue distributions of complex noncentral Wishart matrices. Since in the slow fading scenario it is not feasible to 
determine the ergodic capacity, a metric denominated outage probability is required to evaluate the system performance a. 
The outage probability is related to the cumulative eigenvalue distribution function of a Wishart matrix a, El, nni. Due 
to physical nature of wireless channel and all possible arrangements for antennas arrays, different types of Wishart matrices 
have been studied, such as central and noncentral, associated with Rayleigh and Rician fading, respectively; uncorrelated, 
semi-correlated, and double-correlated, associated with the antenna correlation at the transmitter side and at the receiver side 

El, El- 

B. Extension for MIMO multiuser case 

All the works previously mentioned are concerned with the MIMO single-user channel where the majority of the problems 
are already solved or at least well understood El- However, for the MIMO multiuser scenario there exists many open problems, 
such as the general capacity for the MIMO Relay channel. 

In a wireless multiuser channel, we are generally more concerned in the overall information rate (capacity) of the system 
than the individual user rates as in the single user channel iflTl Ch. 15]. In this way, we could define metrics associated with the 
joint users performance. We have, for example, symmetric capacity and sum capacity. The former is the maximum common 
rate at which both users can simultaneously reliably communicate; the latter is the maximum total throughput that can be 
achieved ||9] pg. 230] and can be seen as a constraint that limits the individual rates of each user. Since sum capacity reflects 
an overall system performance, this metric is of great interest from analytical and practical points of view. 

As can be inferred from the sum capacity name, to evaluate this metric we have to add up the rates of each user. This 
operation leads to a summation of Wishart matrices associated with each one of the MIMO channels involved in the multiuser 
system. Eor example, this situation occurs in two well-known multiuser channels: (i) MIMO multiple access channel (MAC), 
where K users with multiple transmit antennas communicate with one destination also with multiple receiving antennas El; 
and (ii) MIMO Relay channel, where a MIMO transmitter communicates with a MIMO receiver with the help of a MIMO 
relay El- Eor MIMO MAC the sum capacity is a desired metric on performance El- Eor MIMO Relay, the sum capacity is 
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used to determine the cut-set upper bound on the channel capacity IfT^ . 

C. On the Paper Contribution 

Based on our discussion in the preceding section about the sum capacity, hereinafter, we will analyze this metric under 
fast-fading Rayleigh distribution. Therefore, our aim is to determine the ergodic sum capacity for MIMO multiuser scenario. 
Our idea is to use the framework identical to that of a single-user case. It means that we wish to obtain the ergodic sum 
capacity using the marginal eigenvalue distribution of the sum of Wishart matrices. 

For a single-user case, the probability density function of the eigenvalues of a Wishart matrix was given in ©, and since then 
many advances have been achieved for the most variate cases of Wishart distributions. More recently, results on the product 
of rectangular random matrices have appeared in El and El and the authors have investigated ergodic mutual information 
in MIMO communication channel with multifold scattering. By contrast, the progress for the eigenvalues distribution of the 
sum of Wishart matrices has not been going at the same pace. 

Although the most well-known result for the sum of Wishart matrices dates back to 1960’s, it is valid only for the specific 
case where all matrices have the same covariance matrix El; not much is known for the general case of arbitrary covariance 
matrices. In 1201 , the authors have considered linear combination of central Wishart matrices with positive coefficients. They 
have proposed approximating the distributions of the linear combination by central Wishart distributions. Furthermore, in the 
context of multivariate Behrens-Fisher problem, a similar approximation to solve the linear sum of Wishart matrices has been 
given in Ell. Therein the authors have approximated the sum by a single Wishart distribution by determining the associated 
degree of freedom and the parameter matrix. A very recent work in this direction is by one of the present authors, where exact 
matrix distribution has been computed for the sum of two Wishart matrices with arbitrary covariance matrices fT2\ . Moreover, 
explicit result for the eigenvalue statistics has been worked out for the case when one of the Wishart matrices possesses 
covariance matrix proportional to the identity matrix. In the present work we are concerned with the eigenvalue statistics for 
the sum of arbitrary number of central Wishart matrices with covariance matrices proportional to the identity matrix. 

For the MIMO MAC channel, the ergodic sum rate capacity has never been obtained due to the lack of analytical results 
for the joint eigenvalue probability density function of sum of Wishart matrices. Flowever, the capacity with perfect channel 
state information at receiver and transmitter (CSITR) sides is very well studied. With perfect CSIT and CSIR the system can 
be viewed as a set of parallel non interfering MIMO MACs. Thus, the ergodic capacity region can be obtained as an average 
of these parallel MIMO MAC capacity regions (see ifTSll and the references therein). Another approach is to obtain asymptotic 
results on the sum ergodic capacity of MIMO MAC channels. This can be done by considering that the number of receive 
antennas and the number of transmitters tend to infinity El- 

In order to solve these and related challenging problems, we have proposed two distinct approaches that are presented in 
Section [Bll The first approach is the derivation of an exact closed-form expression for the marginal eigenvalue distribution of 
the sum of Wishart matrices. The main idea behind this solution is to demonstrate that the matrix resulting from the weighted 
sum of K Wishart matrices can be rewritten as the product of a single matrix and its conjugate transpose. This resulting 
Wishart matrix happens to correspond to a covariance matrix which incorporates the information about the weights. Therefore, 
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its eigenvalue distribution follows from the pre-existing knowledge about Wishart semicorrelated matrices. The derivation of 


this result is given in Appendices lAl and IbI Our second proposed solution is to approximate the sum of K independent Wishart 
matrices by just one equivalent Wishart matrix. This approach is based on the idea of equating the cumulants, as done in ll20l 
for the case of general covariance matrices. We have found a simple and compact closed-form expression to determine the 
degrees of freedom of this equivalent Wishart matrix. 

In order to show that our proposed solutions are valid, we have chosen the two MIMO multiuser scenario described before, 
viz. MIMO MAC and MIMO Relay. First, by considering an arbitrary set of parameters, we show that the Monte Carlo 
simulated eigenvalue distribution of the sum of Wishart matrices is perfectly described by our exact expression. Then, we 
apply our approximation to find an equivalent Wishart matrix, and compare its eigenvalue distribution with the simulated 
results. The results are promising. 

Moving a step forward, we present in Section |IV] a new closed-form expression for the ergodic sum capacity. This expression 
takes as input the exact eigenvalue distribution of the sum of Wishart matrices or the approximate eigenvalue distribution of 
the equivalent Wishart. We show that the analytical ergodic sum rate capacity matches the simulation results perfectly. All 
these results are shown in Section |V] and give basis for our conclusions presented in Section |Vl] 

Besides all the sections mentioned above, we present some fundamental definitions about Wishart matrices in Section |II] 

II. Preliminaries 

In this section we begin with the definition of complex Wishart distribution, which depends crucially on variance and 
degrees-of-freedom parameters. These are then used to construct the matrix model of our interest, namely the weighted sum of 
central Wishart matrices. This, in turn, is used in later sections for derivation of the probability density function and relevant 
metric for our problem. 


Given a random rui-dimensional non-negative definite matrix with pi degrees of freedom The distribution 

law of Wi, 


Pw,(W,) oc exp (-trS-iW,) , 


( 1 ) 


is called complex central Wishart distribution il, USD, and is denoted by 




( 2 ) 


Here, is the covariance matrix, det(-) and tr(-) represent determinant and trace operators, respectively. In the following we 
will consider = aflrm, where 1^^ is the identity matrix of dimension rrii. 

Consider K independent matrices with the distribution given by (O. We are interested in the eigenvalue statistics of the 
weighted sum of these K matrices normalized by their respective degrees of freedom, viz.. 


K 



( 3 ) 


where Oi G R’*'. Note that the above sum is possible only if the rui’s are identical, say m. 
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It is known for the general case of W = ^ W^. (pi,S), that W ^ Wm* P*; CSl 

Theorem 7.3.2.]. On the other hand, if the covariance matrices S^’s are not proportional to identity matrix, then obtaining the 
distribution of W and its eigenvalues is nontrivial; see for example ll22]| . However, if S^’s are proportional to identity matrix, 
then as shown in appendixW actually corresponds to a semicorrelated Wishart distributed case ll24l . 

For the scenario oc Im, without loss of any generality we may consider = a^lm, as different <7^ values can be 
absorbed in all Let us define 


Vi = {ailpi)a'^, i = (4) 

K 

P = '^Pi- ( 5 ) 

i=l 

With these definitions we now present the exact as well as approximate solution concerning the eigenvalue statistics of W. 

III. Proposed Solution 

This section presents our contribution to determine the eigenvalue distribution for the weighted sum of K Wishart matrices, 
as defined above. First, we present a new exact closed-form expression. Then, we propose an approximation that replaces the 
weighted sum of K Wishart matrices by an equivalent matrix. 


A. Exact closed-form expression for the marginal eigenvalue distribution 


The main result of our paper is given in the following proposition. 

Proposition 1: The marginal density of eigenvalues of W defined in Q is given by 


0 

Pa (A) = cdet 

The entries fj{v,X),gi{\),hij{y) inside the determinant in the above expression are respectively given by 


=1,...,P1 


=l,...,p 

=1,...,P1 

[^iKyjK^'^K)] iK=l 
jK = ly 


( 6 ) 



( 7 ) 

g,iX) = X^-yr{m-t + l), 

( 8 ) 

h- iv) - 

( 9 ) 


Here r(-) is the Gamma function given by r(z) = P * dt, and L^if\x) are the associated Laguerre polynomials ll^ 
eq. 22.5.38]. The normalization c in (|6|l is obtained using 


c 


-1 


—mdet 




i^iKOK {vk)] 


( 10 ) 


Proof: See Appendix A. 

^(T may also be absorbed in a^. Therefore, corresponds essentially to an uncorrelated Wishart case. 
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B. Approximation for Sum of Wishart matrices 

Proposition 2: Given the weighted sum of K Wishart matrices normalized by their respective degrees of freedom as ©, 
we propose the following approximation 


W a; S 




( 11 ) 


where S --- CWm., {ps, = cr^I™, and ps given by 


Ps = 


(eL-.)'' 

2^i—l r)i 


( 12 ) 


with [•] representing the nearest integer operator. 

Proof: The rationality for the approximation is as follow. The expected value of given in (l2]l is given by ||2^ . 


E[Wi] = pi'Sli, 

The variance of the main diagonal elements are given by ll2^ 

var[W,(j,j)] =p^a^, 


( 13 ) 


(14) 


where we have dropped the subscript of a as explained before. 
Also, the expected value of W in © is given by 


E[W] 


K 


Pi 


K 


= Y, ^ = aHm Y ' 


K 


K 


Pi 


and the variance of the main diagonal elements is given by 

K 


var[W(j, j)] = var 


E-w.(j)j) 


K 2 




^ 2 

a: 


i=l 


Pi 


i=l 


Pi 


Notice that the expectation value of RHS of (fTTT i is given by 


E 


K 




and the variance of the main diagonal elements is given by 






4 

= a p, 




= a 


( 15 ) 


( 16 ) 


(17) 


( 18 ) 


Here, we call attention to the similarities of results of (fTSl) and (fTTl l. as well as (fThl) and (fT^ . Therefore, it is possible to 
state a relation between the degrees of freedom of different Wishart distributions by equating (fThl l and (fTSl) and obtaining the 
closed-form expression for ps given by (Ell. Also, since ps is related with the number of columns of the S, it should be an 
integer number, and this is why we have to apply the nearest integer operation in (El- 
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IV. Application 

Generally, a single user communication under fading conditions has a received signal expression given by m (5.86)] 

y — hx + z (19) 

where z ~ CAf {0, 1) is the noise, x is the Gaussian distributed input signal with power constraint ||x|| " < tti, and h is the 
channel gain. Herein, we assume h ^ CAA(0,(T^), therefore, the channel is under Rayleigh fading. Now, suppose that the 
source has Mi transmitting antennas, and the destination has Ni receiving antennas. Hence, the wireless channel is described 
by the complex Ni by Mi random matrix H^, and the received signal is given by S (7.1)1 


y = Hix + z 


( 20 ) 


where z ^ CA/'(0,ImJ is the white Gaussian noise vector at a symbol time (not show here by simplicity), x G CAf^\ and 
y G CN^' . The entries of are h^t, with 1 < r < and 1 < f < Mi. Define a matrix as HI 

f h,h| n,<m, 

w, = <^ (21) 

y H|H, N, > 

where | denotes the transpose conjugated matrix operator. Hence, has real, non-negative eigenvalues Il27ll . The matrix 
is distributed as da with Pi — max(Mi, Ni) and mi — min(Mi, Ni). 

Telatar has shown in his canonical paper i), that the ergodic capacity for the system described by (l20l l is given by 


C = Ew,. 


log2 det (I 


Pi 


W, 


= mi 


log2(l + A)Pa(A) d\. 


( 22 ) 


where P\{X) is the the marginal density of eigenvalues. Since in this work we are interested in multiuser scenario instead of 
a single user described above, we should adapt the capacity equations for a multiuser case. 

The first scenario is the MIMO MAC depicted in Fig. [T] The channel capacity for a MIMO MAC network with K sources 
and one destination under Rayleigh fading is given by m 

/ K 


Cmac = Ewi 


log2 det [ Im + X! 


= E^[log2det (I™+W)] 
= E:^ [tr log 2 (I™ + W)] 


(23) 


where we have used first the equality given in Q, and then the property of matrices that asserts that det (exp (A)) = exp( tr (A)) 

ED. 


In order to solve Cmac given in (|2^ . we can use the eigenvalue distribution given in (|6]l to obtain a closed-form expression 
for the ergodic sum rate capacity given in (l25l) with the appropriated parameters of the Wishart matrices given in the problem 







Fig. 1. System Model of MIMO MAC channel. The destination (D), equipped with receiving antennas, receives signal from K sources, each one 
equipped with Mi transmitting antennas, with i = 1,.. ., K. 


statement. In terms of the marginal density P\{X) of eigenvalues of W, we can write 


Cmac = m Px{X) log2(l + A) dX. 


(24) 


As shown in the appendices, an exact closed form expression for the ergodic capacity can be obtained in determinantal form 
involving Meijer-G functions, and is given by 


Cmac = -m c ^ det 

fl^l 




,pi 


where {v) is given by 




(25) 


with 




(u) 


{ g^,Jiv), i = n, 


Gijiv) = 


, 4-1 


G 


-3,2 


(ln2) r(m — i + 1) 


/ 


V 


0, i — 1, J 1 

1-1 ■ 1 V 

1, z - 1, m, j - 1 


(26) 


(27) 


and hij{v) as in 
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Fig. 2. System Model of MIMO relay channel. The source (S), equipped with Ms transmitting antennas, wants to communicate with the destination (D), 
equipped with receiving antennas. The relay (R), equipped with Mr and Nr transmitting and receiving antennas, respectively, could collaborate in this 
communication. 


On the other hand, by using the proposed approximation given in (fTTl) , (1231) becomes 

/ K 


Cmac = Ewi 


' Ec 


log 2 det ^ 

log2det + 


(28) 


Hence, we can use (l25l l with pi = ps and ai = '^^en set X = 1 to obtain Cmac- 

Now, let’s turn our attention to the MIMO Relay channel shown in Fig. |2] This channel can be viewed as a composition 
of a MAC and broadcast channel (BC) with K = 2. Suppose that the channel gains are known at the corresponding receivers 
only (CSI). In this scenario, an upper bound on the ergodic capacity of the MIMO relay channel is given by HU Theorem 
4.1] 

Clipper — min (Cbc , Cmac ) (29) 


where 


Cbc = Ewi 


log 2 det I, 


Pi 


Wi 


P2 


W; 


(30) 


and Cmac as in (l2Tt . Notice that Cbc is similar to Cmac given in (l2Tt . Therefore a similar procedure can be implemented to 
obtain Cbc- 


V. Numerical Results 

In this section we have obtained numerical results for the closed form expressions and for the proposed approximation. The 
results are compared with Monte Carlo simulations to validate the analytical expressions. For each one of the simulations, 
40,000 channel realizations were performed. In all cases, there is a perfect agreement between analytical and simulation results. 
We have chosen three arbitrary scenarios and one well-known scenario from M- 

Consider a MIMO MAC scenario shown in Fig. [T] with K = 5 users, each one with Mi = 4 transmitting antennas, where 
i = 1,... ,K. Destination node D has Nd = 4 receiving antennas. The normalized signal to noise ratios (ai) at destination 
were arbitrarily chosen, and are shown in Table |I] (Case I). The marginal eigenvalue distribution is shown in Fig. |3] Notice the 
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TABLE I 

Simulation Parameters and Results eor C 


Case 

Oj (dB) 

MIMO 

Ps 

Ergodic Sum Rate Capacity (bits) 

Simulation 

Analytical 

Approximation 

I 

ai = 19.8 
a 2 = 29.5 
a 3 = 29.8 
a 4 = 26.1 
= 21.7 

4x4 

13 

44.20 

44.20 

44.15 

II 

= 28.3 
02 = 17.7 
03 = 26.5 
04 = 27.3 
05 = 29.3 
06 = 21.5 
07 = 19.5 
og = 27.9 
Og = 9.3 
oio = 24.3 

8 x 8 

51 

93.86 

93.86 

93.85 

III 

oi = 9.7 

02 = 17.6 
03 = 16.7 

2 x 2 

3-4 

10.94 

10.94 

11.01 



Fig. 3. Marginal distribution of eigenvalues of W for Case I (see Table H). Simulation results are in perfect agreement with the closed form analytical 
distribution (red line). The distribution for the proposed approximation, shown with blue dashed line, is also very close to the exact result. 


perfect agreement between the simulation and analytical results. The approximate ergodic sum rate capacity was also computed 
using the equivalent matrix S (with ps = 13 degrees of freedom) and is shown in the last column of Table U The eigenvalue 
distribution for the approximation is also shown in Fig. [3 Note how close the approximation and the ergodic sum rate capacity 
are to the exact values. 

In the next scenario we have increased the number of users to if = 10 with MIMO 8 x 8 . The coefficients are given 
in Case II of Table U Notice again in Fig. |4] a perfect match between analytical and simulation results. The ergodic capacity 
results also agree, as can be seen in Table |I] The eigenvalue distribution from the approximation is shown in Fig. |4] as well. 

Now let us move on to a more involved scenario depicted in Fig. |2 The ergodic capacity was originally calculated in 
IfT^ using convex programming. The parameters used are depicted in Case III of Table |T] and the plot for the eigenvalue 
distribution is given in Fig. 0 Fig. | 6 ] shows the eigenvalue distribution for the MAC channel with the following parameters 
(o 2 = 17.6, as = 16.7). The upper bound on ergodic capacity, as mentioned before, is the minimum of the capacity of BC 
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A 


Fig. 4. Marginal distribution of eigenvalues of W for Case II (see Table Q). Simulation results match perfectly with the closed form analytical distribution 
(red line). The distribution for the proposed approximation is also shown using blue line. 



Fig. 5. Marginal distribution of eigenvalues of W for Case III (BC). Simulation results are in perfect match with the closed form analytical distribution 
(line). 


and MAC, as given in Table U 

Besides the Case III, we have also reproduced the scenario given in ifTbl Fig. 5] for MIMO Relay channel. This scenario is 
well known because the upper bound and the lower bound “converge”. That is to say, the ergodic capacity of the MIMO relay 
channel over Rayleigh fading can be characterized under this SNR condition HE). The constraints for this scenario are 02 = 03 , 
oi = 1002 , and 0 < 02 < 30 dB. The ergodic capacity results for Monte Carlo simulation and the proposed approximation 


0.008 

0.006 

cC 

0.004 

0.002 

0.000 

0 50 too 150 200 250 300 350 

A 

Fig. 6. Marginal distribution of eigenvalues of W for Case III (MAC). Simulation results match perfectly with the closed form analytical distribution (line). 
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Fig. 7. Ergodic Capacity Upper Bound versus a 2 with constraints as = a 2 and ai = 10a2. Based on Ql Fig. 5]. The simulation results are in perfect 
agreement with the analytical results obtained with the proposed approximation. 



Fig. 8. Error of approximation versus weight ratio ai/a 2 and a 2 = 5. 

are given in Fig. |2l Notice again the perfect agreement of the results. 

The results in Table U and in Fig. [T] show that the proposed approximation results in a very small difference from the exact 
result for all mentioned scenarios. Since the approximation depends on the weight factor of each matrix, or in other words 
depends on the SNR of each channel, it would be interesting to investigate cases where these weight factors varies. Fig. 0 
shows the percentage error as the ratio 01/02 varies from 0 to 15 dB. Note that the error is less than 1% for 0 < 01/02 < 2 
and 01/02 > 7. 


VI. Conclusions 

In this work we uncovered a one to one correspondence between the weighted sum of arbitrary number of uncorrelated central 
Wishart matrices and a single semicorrelated Wishart matrix. Using this observation we presented a closed form expression for 
the marginal distribution of the eigenvalues for the weighted sum of K complex central Wishart matrices. To the best of our 
knowledge this problem has not been tackled before. Here the motivation for establishing a result emerged from the multiuser 
information theory area. However, since Wishart matrices play crucial role in diverse helds, we believe that our results are 
relevant to these as well. We would like to remark that it is also possible to obtain results for the joint probability density of 
all eigenvalues, and correlation functions involving distribution of two or more eigenvalues. 

We applied our new closed-form expression for analyzing the ergodic sum rate capacity of MIMO multiuser channels. 
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Moreover, we also derived a closed form expression for the ergodic channel capacity and used it to obtain the capacities for 
MIMO MAC and MIMO Relay channel. Besides the closed form exact expression for marginal distribution, we also proposed 
an approximation that is very simple and presents promising results when used to obtain ergodic sum capacity. In addition, 
we confirmed the validity of all our analytical expressions by Monte Carlo simulations. 


Appendix A 

Mapping to a semicorrelated Wishart distribution 
Consider pi x m dimensional complex matrices H^, i = 1,..., AT, from the normal distribution: 

Pi(Hj) oc exp 


trH,Sr'Hn , 


(31) 


where Then m x m-dimensional matrices = H|Hi, are respectively from the complex-Wishart distribution 

given in (|2]i. The matrix W can be written as 

2=1 2=1 


Hi nl H 


K 


—I 0 

Pl^ 


0 


P2 ^ 


= tfDH = G^G. 


O-K f 

Pk 


■Pk 


Hi 

H2 


H 


K 


(32) 


We defined here 


Hi Hi 


Ht = 

D = diag 

iPi 

G = 


ai j 


■Pl 5 


Hi, 


) ^PK 

PK 


(33) 

(34) 

(35) 


With the above information it is clear that G satisfies the distribution 


Pg(G) oc exp[-tr(GtV-iG)], 


(36) 


where V = diag [uilpj,..., with Vi (and p) as defined in (|5]l. Therefore, we are looking essentially at a semicorrelated 

Wishart case, with the diagonal-covariance matrix possessing some equal-value entries (multipUcities/degeneracies). Thus, the 
problem boils down to determining the eigenvalue statistics of G^G. We start with the case of a diagonal covariance matrix 
with unequal entries along the diagonal, i.e., we will use V = diag (vi,..., Dp) in the above distribution instead of V, work 
out the results for this case, and eventually take adequate limits to obtain the case of V. 

For the semicorrelated Wishart matrices, exact result for the marginal density of eigenvalues is available from several notable 
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works Il29l - ll^ . We use here the form derived in OTlI . |[32l. Consider p x m dimensional complex matrices G taken from 
the distribution (l36l l. but with the covariance matrix V as defined above. The marginal density of p eigenvalues of GG^ for 
p < TO, is given to be 

exp( —A/'O. 


PxW = - 


1 


pAp({i)-i}) 


det 


0 




vr ■ 1 


(37) 


[r(m—z+ 1 ) 

Here Ap({{i“^}) = ~ Vandermonde determinant. Since G^G and GG^ share the same 

nonzero eigenvalues, the above result holds for p > to as well, with the factor p appearing in the denominator of the prefactor 
replaced by to. Moreover, in this case the bottom p — to entries of the first column in the determinant comprise diverging 
gamma function in the denominator, and hence become zero. 

We now use the relation (l22l l to derive the ergodic channel capacity. To this end we expand dTTl) using the first column and 
obtain 

1 


Pa(A) = - 


FTTT E(-l) 


^ -det 


H) r(TO-p +1) 


exp( —A/vj) 






(38) 


Next we bring in the A™ ^/r(TO — p + 1) factors occurring before the determinants to the respective first rows, i.e., with 
exp{—X/vj)/Vj^. This gives 


^a(A) = - 


1 


3TTT 


toAp({i)-i}) 


^m-f^ exp( — A/vj 


r(m-/j+l) 




\y^ 


(39) 


The above equation serves as yet another expression for the marginal density. 

For ergodic capacity we use (l22l l and obtain the following expression by interchanging the A-integral and the summation; 


/ 


C = — 


iPnTE(-i)‘ 


d\ 


det 


The A-integral can be brought into the first row of the determinant, along with the factor log 2 (l + A) to yield 


V 


— exp( —A/{) 


r(m—/T+ 1 ) vV 




[Vj *'''^] = 




log2(l -f A). 


(40) 


C = - 






[Vj = 


where 


Gp,3iv) = 


A™-'" exp(-A/%) 


dX 


r{m — P + 1) v)p 


log2(l + A). 


(41) 


(42) 


This integral can be expressed in a closed form with the aid of Meijer-G function If33ll . This is facilitated by considering the 
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following special cases of Meijer-G functions: 

/ 


G, 


1,0 

0,1 


/3 


z = z^e 


G. 


2.2 


1 , 1 
1,0 


z = ln(l + z), 


(43) 


We also use the convolution integral satisfied by Meijer-G function: 



( 

\ 


( 

\ 

poo 

Ol, • • • dp 



Cl, ' ■ ' C(j 

1 

/ dzG^;P 


pz 



UJZ 

h 

^ 61, ••• 6, 

J 


y dl, ■■■ dr 

) 


_ 1 ^n+fi,7n-\-u 

~ ^q+<T,p+T 


= -G: 


p+T,q+a 


til, * * * , Uyj, di, * * • d-T ^ iiyj_|_i, * * * , U. 
Ul, * * • , dn , dl, ■ • ■ dr , Un-t-l, ‘ ‘ , dp 

5 * * ’ 7 Cl j * * * Ccr, 6772-1-1, * * ' 7 6g 



(44) 


The restrictions on the indices for this integration formula can be found in 1^ . Therefore, we obtain a closed form expression 
for Qij{v) as given below in (l47l i. Afterwards we perform row interchanges in the determinants to bring Qpj in the respective 
/ith row. This leads to the removal of (—I)'' factor. Consequently, we arrive at the following expression for ergodic channel 
capacity: 


Here 


C = -- 


det[' 0 . 




E 


det 




J i,j=l,...,p 


(45) 


with 








1 

(In 2) r(r7i 


i + 1) 


*-^ 2,3 


( ■ 1 ■ 

2 — 1,7 

1) 

^ 2 — 1, 2 — 1, m 

Vk 1 


(46) 


(47) 


Appendix B 

Proofs for equations dUi and (|25T i 

To obtain equations (|6]) and (|25]) we need to assign vi = ... = Vp^ = vi]Vp^j^\ = ... = Vp^ = V 2 ;■■■ ; '0p(7^_i,-|-i = ... = 
'^PK — fA ill (tT7l) and (l45l) . However, direct substitution of these values makes the determinant in the numerator, as well as 
the determinant in the denominator to become zero. Therefore, we must carry out this substitution in a limiting manner, as 


described below. 




















16 


Let us pay attention to the columns involving up to pi in ( iJTl i. The ratio of the determinants appears 


det 

0 

( 0 -l)mg 2 )/A (^-l)mgi) 2 ^A (t)+)™e 

... ( 1 )+)- 

D ^ A 

PI 

-1 



r(m— 2 + 1 ) 



det 

[v-y-^ (u-i)-i ... (++-1 ... 



(48) 


Consider for j = 2, v, = ^ with small 6j, and Taylor-expand up to <5^ ; 


{vj 


1 ^ ^ 

lx Cl 1 


r—0 




r=0 


r! 


Now, applying adequate column operations we obtain 


det 




_ 

IFT)! a{vYY 




r(m— 2 + 1 ) 




9(«i 




<5^- 


(/-!)! d{vYy 




det 


•'^1 > ) •■• 0-1)! a(ifTj-rTlt'l ) 


sr 


— l\z—1 


(49) 


The factors containing 6j and factorial can be canceled out after being pulled out of the columns, both from numerator and 
denominator. Therefore, we are left with 


det 


0 


r(m—2+1) 




QA-iy-i 


det 


(^1 ) a(v-^)^^i ' ■■■ ' 


(50) 


We have +_i i (0^ ^)® ^ = (r(i)/r(i — r + l))('Ci The derivatives of (v - can be evaluated with the aid 

0(Vj^ 

of Rodrigues’ formula for the associated Laguerre polynomials, 

, _z-V 5^= 


irw = 


k! dz^ 




( 51 ) 


using adequate scaling of the variables. By implementing similar steps for rest of the columns, we arrive at (HJ. 
Similar steps can be used to arrive at ( l25l l, starting from ( l45l l. The derivative of Meijer-G follows from the result 


dz^ 


/ 

\ 


f n 

\ 

ai, 02 


= g^;4 

0 , Ol, 02 


&1+2+3 



^ bi, & 2 , bs, r 

) 


which is a special case of the following more general identity 

/ 


dr 

r ^m,n 

dz^ p’’^ 


Q-l, • • • , Qjp 
bi, ■■■ ,bg 


_ ^ 771 , 71+1 

— ^p+l,(?+l 


0, oi, ■ • ■ , ap 

bi, ■■■ , b„, r 


(52) 


( 53 ) 
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